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Recent thermodynamic measurements on two-dimensional (2D) electron systems have found di- 
verging behavior in the magnetic susceptibility and appearance of ferromagnetism with decreasing 
electron density. The critical densities for these phenomena coincide with the metal-insulator transi- 
tion recorded in transport measurements. Based on density functional calculations within the local 
spin-density approximation, we have investigated the compressibility and magnetic susceptibility 
of a 2D electron gas in the presence of remote impurities. A correlation between the minimum in 
the inverse capacitance {dfi/dn) and the maximum of magnetization and magnetic susceptibility is 
found. Based on values we obtain for the inverse participation ratio, this seems to be also the MIT 
point. 

PACS numbers: 71.10.-w, 71.10.Ca, 71.30.-|-h 



Introduction-There has been a large amount of exper- 
imental and theoretical activity in recent years to un- 
derstand the ground state properties of two-dimensional 
(2D) electron systems^. Most notably, the observation of 
a metal-insulator transition (MIT) in these systems pro- 
vides a major motivation to study the various physical 
properties. Most experiments perform transport mea- 
surements obtaining resistivity or conductivity as a func- 
tion of temperature at varying electron density to de- 
duce the metallic or insulating phases. In contrast, Hani 
et al^ used the capacitance technique, a thermodynamic 
measurement, to measure the compressibility finding that 
at high densities the results for the inverse capacitance 
{dfi/dn) agree with Hartree-Fock or mean- field theories, 
whereas below the critical density the system becomes 
very inhomogeneous and dn/dn goes up as the density is 
lowered. Similar observations were also made by Dultz 
and Jiang'^ on a 2D hole system who noted that the in- 
verse compressibility is minimum at the same density 
where the MIT occurs. 

Recent experiments at in-plane magnetic field concen- 
trated on the spin susceptibility, Lande ij-factor, and ef- 
fective mass of the 2D electron systems present in Si- 
MOSFETS and GaAs quantum-well structurea'^i^^i'^i^i^. 
In particular, Shashkin et aliSilL reported a sharp in- 
crease of the effective mass near the critical density at 
which the system starts to show deviations from the 
metallic behavior. On the other hand, Pudalov et alA 
have found only moderate enhancement of the spin sus- 
ceptibility and effective mass in their samples. Ther- 
modynamic measurements of magnetization of a dilute 
2D electron system were reported by Prus et al.^ and 
Shashkin et ai?- Both experiments found large enhance- 
ment of the spin susceptibility Xs over its Pauli value. 
Whereas the measurements of Prus et alfi found no in- 



dication toward a ferromagnetic instability, Shashkin et 
al.^ observed diverging behavior in Xs at a critical density 
coinciding with the MIT density obtained from transport 
measurements. 

On the theoretical side, calculation of compressibil- 
ity for a 2D system of electrons in the presence of 
disorder predicts the observed behavior of upturn and 
divergencelS.. Shi and Xiei2. performed density functional 
calculations based on the (unpolarized) local density ap- 
proximation developed by Tanatar and Ceperlej*i^ within 
the Thomas-Fermi (TF) theory, and found similar results 
for the compressibility. They also identified the MIT 
point with the percolation transition point in this system. 
This theory was further strengthened by Das Sarma et 
al.^^ who measured the critical exponent for the conduc- 
tivity and found it in agreement with that proposed by 
percolation theory. 

In this work, we investigate the spatial distribution of 
carrier density and magnetization, as well as the com- 
pressibility and spin susceptibility of a 2D electron sys- 
tem using the local spin-density approximation (LSDA) 
both at the TF and Kohn-Sham (KS) levels. Localiza- 
tion properties are characterized by the inverse partic- 
ipation ratio (IPR). The functional we use is the one 
constructed by Attacallite et alt^ from the very recent 
quantum Monte Carlo calculations for correlation energy 
appropriate for uniform systems. An important feature 
of these simulations is that a transition to a ferromag- 
netic phase at low densities is built in the functional. A 
disorder potential due to remote impurities is included to 
make the calculation realistic. The density distribution 
of the system shows high and low density regions as re- 
ported earlier by Shi and Xie^'^ which supports the idea of 
percolation transition. We also find correlated behavior 
of inverse capacitance minimum and spin susceptibility 
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maximum, which gives thermodynamic evidence for the 
appearence of a phase transition. Thus, a unified picture 
for the thermodynamic behavior of a dilute system of 2D 
electrons in the presence of disorder emerges from our 
calculations. 

Theory-We consider a 2D electron system interacting 
via the long range Coulomb interaction Vq — 2Tre^/{eq) 
where e is the background dielectric constant. The sys- 
tem is characterized by the dimensionless interaction 
strength = l/^nnag'^y^'^, where n is the 2D electron 
density and a*^ = ?i^e/(m*e^) is the effective Bohr radius 
defined in terms of the band mass m* of electrons in the 
bulk semiconductor structure. 

Within the spin-density functional theory the total en- 
ergy of an TV-electron interacting system in a local ex- 
ternal potential Vcxtii^) is a unique functional of spin- 
dependent densities 7T,|(r) and ni{r). The total energy 
functional can be expressed as 

^^Totai[»^T'"i] = ET[n^,ni] + EH[ni,ni] 

We approximate the kinetic energy functional by the 
Thomas-Fermi- Weizacker (TFW) form given bj*ii 



ET[nT,ni] = ^ J dr 



1 |Vn,(r)|2 



The direct Coulomb energy is given by 

/■ 9 9 , n(v) n(r') 
EH[n^,n,] = - J d2rdVAi_A_i 



(1) 



(2) 



where n(r) = n|(r) + n|(r) is the total density. The 
local spin-density approximation proposed by Attaccalite 
et al.^ is used to calculate the exchange and correlation 
potentials. 

The disorder studied in this work comes from a random 
distribution of charged impurities with density rii and at 
a setback distance d from the electron layer. The energy 
functional due to the external potential is 



Ecxt [n] 



d^rl/oxt(r) n(r) 



(3) 



where the external potential, due to remote impurities, 
is given by 



14xt(r) = -^ 



Ze^e 



[(r-ri)2 -f-d2]i/2 



(4) 



Method- The spin-densities nf(r) and nx(r) that ex- 
tremize the total energy functional can be obtained by 
annealing from a Monte Carlo (MC) simulation. A suf- 
ficiently high temperature is first chosen and a Metropo- 
lis MC run is performed long enough to reach thermo- 
dynamic equilibrium. Then the temperature is reduced 
and the run is repeated. This is continued until the final 
temperature is sufficiently low so that very little energy 



fluctuations occur during the last run. This simulation 
is done in order to reach the global minimum of the en- 
ergy landscape. Once one is near the bottom of the val- 
ley, eventually a steepest-descent algorithm is applied to 
reach the minimum energy structure faster. The areal 
integral in the energy functional is approximated by a 
discrete sum: basically the density and potentials are 
discretized on a 44 by 44 mesh whose size is one effec- 
tive Bohr radius being equal to 100 A in GaAs samples. 
The long-range Coulomb potential is calculated using the 
Ewald sum method. 

The results from the TFW kinetic energy functional 
were also compared to the solution of the KS equations 
where the kinetic energy is calculated exactly. We note 
that a similar Hartree-Fock calculation was recently per- 
formed on a Mott insulator in the half-filled limit^^ where 
the ground state is antiferromagnetic. We believe that 
the exchange-correlation potential we are using is more 
realistic, and contains the correct physics, namely that 
at very low densities, before the Wigner crystallization, 
the clean system becomes ferromagnetic. 

Results and Discussion-For a given distribution of the 
impurities, the total energy is computed as a function of 
density. The obtained curve was then fitted with 3 best 
analytical curves possible and then derivatives were taken 
analytically in order to reduce fluctuations. The inverse 
compressibility is obtained as = v?d^/dn. There is 
some variations in the results at higher densities, but as 
Fig. 1 shows, both the minimum and the zeros of dpi/dn 
are identical with reasonably good accuracy. 



(a.u.) 




Density (10 ^ a. u . ) 



FIG. 1: (Color online) The derivative of the chemical poten- 
tial versus density, expressed in units of the effective Bohr 
radius, obtained with the TFW functional. The disorder po- 
tential corresponds to a spacer layer of d = 3 a.u. and a 2D 
impurity concentration of Ui — 2.6 x 10^'^ a.u. both expressed 
in units of the effective Bohr radius. 



The electron gas was found to have a spin polariza- 
tion f(r) = [n|(r) — ni{r)]ln{r) reaching a maximum 
of about 0.05 at B = in the low density regions. 
At lower densities below the minimum in the curve for 
dii/dn the electron gas separates into electron puddles 
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where the disorder potential is lowest. Above that crit- 
ical density there seems to be percolation as noted by 
Shi and Xie4^ The magnetization per electron (^(r)), 
as can be observed in Fig. 2, was found to be largest 
at the critical point (minimum of d^/dn) and it was 
found that both ^ and the susceptibility are highest in 
the low density regions where the depth of the effective 
potentia^Kff = Kxt + ^Hartree + 14c) IS uiost shallow. 




FIG. 2: (Color online) The density, magnetization per par- 
ticle ^, and effective potential distribution at the critical 
point (minimum of d^/dn), within the TFW theory. The 
parameters for this figure are as follows, d = 3a.u., Ui — 
2.6 X 10~^ a.u., rs = 8, n = 5 X 10~^ a.u. One can see that 
the magnetization is largest in the low density or saddle point 
(for the potential) regions. 

To investigate further the effect of the disorder on the 
spontaneous magnetization of the electron gas, we calcu- 
lated its magnetic susceptibility by applying an in plane 
magnetic field to the sample and measure x — 
at different densities. The effect of the parallel field is 
only to create a Zeeman splitting in the energy of elec- 
trons but the field does not affect their orbital motion as 



the system is purely two-dimensional. As can be seen in 
Fig. 3, even for low fields, both the magnetization and the 
susceptibility reach their maximum at the density where 
dji/dn is minimum. 

B B^Omi Ht.ii.i 

' B^0.002(aM.! 

^ " B-0.003 (tLu.) 

. _ I* B=Q.0O5 (a.u.) 

s B=0.007 (tLU.) 




2 4 6 8 10 12 

Deiisity(10- a.u.) 




FIG. 3: (Color online) The spin polarization ^ (a) and its 
derivative x (b) versus the density for different strengths of 
the applied parallel field. Impurity parameters are the same 
as in previous figures. It is seen that the maximum in the 
magnetization or x occurs exactly at the density where d^/dn 
was minimum. 

At its minimum, we have d^/dn < 0, and according 
|-g3ji3 g^j-g ^j^g percolatiou transition. An added elec- 
tron will mostly extend in the percolation region where 
the density is very low and effective potential fiat, and 
according to our exchange-correlation (XC) functional at 
low densities, the added electron tends to become magne- 
tized. This tendency is strongest only at the percolation 
point since at high densities the XC effects do not fa- 
vor magnetization, and below percolation densities the 
added electron will tend to be localized in the puddles 
where again the density is high. To confirm this picture, 
we need to know the density distribution of the high- 
est occupied and lowest unoccupied states, below and 
above the Fermi level respectively, where the subtracted 
or added electrons localize. TFW approximation can not 
solve for energy levels and eigenstates, but one can ex- 
tract these levels from the solutions of the KS equations. 
For the same disorder parameters as above, but for a 
twice coarser mesh, we have solved the KS equations and 
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we are displaying in Fig. 4 the highest occupied (HOS) 
and lowest unoccupied (LUS) states at the critical den- 
sity where B^pildr? — 0. 




FIG. 4: (Color online) The density distribution of the highest 
occupied, HOS (a) and lowest unoccupied, LUS (b) states just 
below and above the Fermi level, respectively, at the critical 
density where d^/dn is minimum. Fig. 5 shows that they are 
both down-spin states. 



constant value of 3, and below this energy it starts going 
up implying that the HOS seems localized and the LUS 
is delocalized. This confirms that at the chosen density 
where d^/dn is minimum, occupied states are localized 
and unoccupied ones delocalized, meaning that one is 
also at the MIT point. 

Summary~h\ summary, we have investigated the 
ground state density and magnetization distribution of 
an interacting electron system in 2D in the presence of re- 
mote impurities. Within the purely thermodynamic cal- 
culations based on LSDA, we found that the percolation 
transition density where the inverse capacitance d^ijdn 
exhibits a minimum, coincides with the MIT based on 
the IPR results. At the same critical density the mag- 
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FIG. 5: (Color online) IPR versus the energy in a.u. for the 
states around the Fermi level (which is at -0.63 a.u.), at the 
critical density where d^/dn is minimum. 



To see whether the HOS and LUS are localized or not, 
one can compute the inverse participation ratio defined 
as, IPR = / \il)\'^<Pr/ J \^p\^(Pr. For an exponentially lo- 
calized state in 2D, it gives the square of the inverse lo- 
calization length. From the KS equations, we have deter- 
mined the band structure of the finite sample and for the 
considered system, it turns out that both the HOS and 
LUS are in the down-spin state, consistent with a large 
magnetic susceptibility. The IPR is plotted as a function 
of the energy of the states and displayed in Fig. 5. One 
notices that above the Fermi level, it reaches an almost 



netization and magnetic susceptibility display a maxi- 
mum where the system goes into a partially polarized 
ferromagnetic state, and maximum in magnetization are 
correlated. These calculations suggest that the percola- 
tion transition, the ferromagnetic transition and the MIT 
points are all correlated. 
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